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Valuation of Mortality Risk via the Instantaneous Sharpe Ratio: 
Applications to Life Annuities 

Abstract: We develop a theory for valuing non-diversifiable mortality risk in an incomplete mar- 
ket. We do this by assuming that the company issuing a mortality-contingent claim requires 
compensation for this risk in the form of a pre-specified instantaneous Sharpe ratio. We apply 
our method to value life annuities. One result of our paper is that the value of the life annuity 
is identical to the upper good deal bound of Cochrane and Saa-Requcjo (2000) and of Bjork and 
Slinko (2006) applied to our setting. A second result of our paper is that the value per contract 
solves a linear partial differential equation as the number of contracts approaches infinity. One can 
represent the limiting value as an expectation with respect to an equivalent martingale measure 
(as in Blanchet-Scalliet, El Karoui, and Martellini (2005)), and from this representation, one can 
interpret the instantaneous Sharpe ratio as an annuity market's price of mortality risk. 

Keywords: Stochastic mortality; pricing; annuities; Sharpe ratio; non-linear partial differential 
equations; market price of risk; equivalent martingale measures. 

JEL Classification: G13; G22; C60. 
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1. Introduction 

A basic textbook assumption on pricing life insurance and annuity contracts is that mortality 
risk is completely diversifiable and, therefore, should not be priced by capital markets in economic 
equilibrium. Indeed, under the traditional insurance pricing paradigm, the law of large numbers is 
invoked to argue that the standard deviation per policy vanishes in the limit and in practice a large 
enough portfolio is sufficient to eliminate mortality risk from the valuation equation. Although this 
only applies to a large number of claims that are independent random variables, it is commonly 
accepted that, in fact, typical real- world portfolios of life insurance and pension annuities satisfy 
the requisite assumptions. 

Thus, for example, if actuaries estimate that 50% of 65-year-olds in the year 2007 will live 
to see their 85th birthday in the year 2027, then a sufficiently homogeneous group of 65-year-olds 
should be charged 50 cents per present value of $1 dollar of a pure endowment contract (ignoring 
investment income). The underlying assumption is the existence of a well-defined population 
survival curve, which determines the fraction of the population living to any given time. This is 
also the probability of survival for any given individual within the group. Thus, if the insurance 
company sells enough of these policies and charges each policyholder the discounted value, on 
average the company will have enough to pay the $1 to each of the survivors. Using the language 
of modern portfolio theory, the idiosyncractic risk, as measured by the standard deviation per 
policy, will go to zero if insurers sell enough policies, so the mortality risk is not priced. 

However, a number of recent papers in the insurance and actuarial literature challenge this 
traditional approach within the framework of financial economics. This research falls under the 
general title of valuation of stochastic mortality. In this paper, we provide a unique contribution to 
this literature by valuing stochastic mortality using techniques that have been traditionally applied 
to portfolio management. 



2 



Within the topic of stochastic mortahty, recent papers by Milevsky and Promislow (2001), 
Dahl (2004), DiLoernzo and Sibillo (2003), and Cairns, Blake, and Dowd (2006) proposed specific 
models for the evolution of the hazard rate. Others, such as Blake and Burrows (2001), Biffis 
and Millossovich (2006), Boyle and Hardy (2004), Cox and Wang (2006), proposed and analyzed 
mortality-linked instruments. Few, if any, have focused on the actual equilibrium compensation for 
this risk. Milevsky, Promislow, and Young (2006) provided some simple discrctc-time examples, 
while Dcnuit and Dhaene (2007) used comonotonic methods to analyze this risk. A practitioner- 
oriented paper by Smith, Moran, and Walczak (2003) used financial techniques to justify the 
valuing of mortality risk, although their approach is quite different from ours. 

What this literature essentially argues is that uncertainty regarding the evolution of the instan- 
taneous force of mortality induces a mortality dependence that cannot be completely diversified by 
selling more contracts. This phenomenon induces a mortality risk premium that is valued by the 
market and whose magnitude will depend on a representative investor's risk aversion or demanded 
compensation for risk. In other words, if there is a positive probability that science will find a 
cure for cancer during the next thirty years, aggregate mortality patterns will change. Insurance 
companies must charge for this risk since it cannot be diversified. The question we confront is: 
how much should they charge? 

From a technical point of view, in this paper, we value mortality-contingent claims by assuming 
that the insurance company issuing the claims will be compensated for aggregate mortality risk 
via the instantaneous Sharpe ratio of a suitably defined hedging portfolio. Specifically, we assume 
that the insurance company chooses a hedging strategy that minimizes the local variance, as in 
Schweizer (2001a), then picks a target ratio a of net expected return to standard deviation, and 
finally determines the corresponding value for any given mortality-contingent claim that leads to 
this pre-determined a. See Bayraktar and Young (2008) for further discussion of this methodology 
from the standpoint of the mean-variance efficient frontier. 

The main contribution of this paper is to develop a financial theory of valuing mortality- 
dependent contracts under stochastic hazard rates. We apply our method to value life annuities. 
Arguably, the primary result of the paper is that the value per contract solves a linear partial dif- 
ferential equation as the number of contracts approaches infinity. One can represent the limiting 
value as an expectation with respect to an equivalent martingale measure, and from this represen- 
tation, one can interpret the instantaneous Sharpe ratio as the annuity market's price of mortality 
risk. 

Although our valuation mechanism is very different-we use a continuous version of the stan- 
dard deviation premium principle-it turns out that the value of the annuity corresponds to the 
upper good deal bound of Cochrane and Saa-Requejo (2000) and of Bjork and Slinko (2006). The 
no-arbitrage price interval of a contingent claim is a wide interval, and Cochrane and Saa-Requejo 
(2000) find a reasonably small interval for prices by ruling out so-called good deals by putting a 
bound on the absolute value of the market price of risk of the Radon-Nikodym derivative of valu- 
ation measures with respect to the physical measure. The lower and upper prices for the resulting 
interval are obtained through solving stochastic control problems, although this is not precisely 
stated in Cochrane and Saa-Requejo (2000). Our valuation method provides a new interpretation 
of the upper good deal bound of a life annuity as the value to a seller under the instantaneous 
Sharpe ratio a when the hedging strategy is chosen to minimize the local variance. The lower good 
deal bound is the value to a buyer of a life annuity. 

The remainder of this paper is organized as follows. In Section 2, we present our financial 
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market, describe how to use the instantaneous Sharpc ratio to value the hfe annuity, and derive 
the resulting partial differential equation (PDE) that the value solves. We also present the PDE 
for the value a^"^^ of n conditionally independent and identically distributed life annuity risks. 

In Section 3, we compare our valuation method with two that are common in the literature. 
Specifically, (1) we consider the good deal bounds of Cochrane and Saa-Requejo (2000) and of 
Bjork and Slinko (2006). We show that for the issuer of a life annuity, our value is identical to the 
upper good deal bound; for the buyer of a life annuity, our value is identical to the lower good deal 
bound. Therefore, we give an alternative derivation of the good deal bounds for our setting. Then, 
(2) we consider indifference pricing via expected utility, as described in Zariphopoulou (2001), for 
example. We observe that the relationship between our valuation method and indifference pricing 
is similar to the relationship between the standard deviation principle and variance principle in 
insurance pricing (Bowers et al., 1986). 

In Section 4, we present several properties of a^"'^ and find the limiting value of ^a*-"-*- We 
show that this limiting value solves a linear PDE and can be represented as an expectation with 
respect to an equivalent martingale measure. Section 5 concludes the paper. 

2. Instantaneous Sharpe Ratio 

In this section, we describe a life annuity and present the financial market in which the issuer of 
this contract invests. We obtain the hedging strategy for the issuer of the life annuity. We describe 
how to use the instantaneous Sharpe ratio to value the life annuity and derive the resulting PDE 
that the value solves. We also present the PDE for the value a^"^ of n conditionally independent 
and identically distributed life annuity risks. 

2. 1 Mortality Model and Financial Market 

We use the stochastic model of mortality of Milevsky, Promislow, and Young (2005). We 
assume that the hazard rate Aj (or force of mortality) of an individual follows a diffusion process 
such that if the process begins at Aq > A for some nonnegative constant A, then \t > \ for all 
t > 0. From a modeling standpoint, A could represent the lowest attainable hazard rate remaining 
after all causes of death such as accidents and homicide have been eliminated; see, for example, 
Gavrilov and Gavrilova (1991) and Olshansky, Carnes, and Cassel (1990). 

Specifically, we assume that 



in which {W^^} is a standard Brownian motion on a filtered probability space (fi,.?^, (.7^t)t>o, P). 
The volatility a is either identically zero, or it is a continuous function of time t bounded below by 
a positive constant k on [0, T]. The drift /i is such that there exists e > such that if < A — A < e, 
then /x(A, t) > for all t G [0,T]. This condition on the drift ensures that At > A with probability 
1 for all t G [0, T\. In general, we assume that jj, satisfies requirements such that (2.1) has a unique 
strong solution; see, for example, Karatzas and Shreve (1991, Section 5.2). Note that if cr = 0, 
then At is deterministic. 

Example 2.1 Suppose /x(A,t) = g + \cr^ + mln(Ao — A) + mgt — mln(A — A) and a{t) = cr in 
which 5 > 0, m > 0, and cr > are constants. The solution to (2.1) is given by 



d\t = n{Xt, t) (At -X)dt + a{t) (At - A) dW/, Aq > A 



(2.1) 



At = A + (Ao - A) exp {gt + Yt) , 
dYt = -mYtdt + a dW^^, Yo = 0, 



(2.2) 
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a mean-reverting Brownian Makeham's law (Milevsky and Promislow, 2001). Note that Yt = 
a e"™*^*"*) dWg, which reduces to aW^ when m = 0. As an aside, /x(A, t) satisfies the conditions 
in Assumption 4.1 below. □ 

Suppose an insurer issues a life annuity to an individual that pays at a continuous rate of $1 
per year until the individual dies or until time T, whichever occurs first. Throughout this paper, 
we assume that the horizon T is fixed. In Section 2.2.2, to determine the value of the life annuity, 
we will create a portfolio composed of (1) the obligation to pay this life annuity, of (2) default-free, 
zero-coupon bonds that pay $1 at time T, regardless of the state of the individual, and of (3) money 
invested in a money market account earning at the short rate. Therefore, we require a model for 
bond prices, and we use a model based on the short rate and the bond market's price of risk. 

The dynamics of the short rate r, which is the rate at which the money market increases, are 
given by 

drt = birt,t)dt + cirt,t)dWt, (2.3) 

in which b and c > are deterministic functions of the short rate and time, and {Wt} is a standard 
Brownian motion with respect to the probability space {0,,J-, (.7-f)f>o, P), independent of {1^/*}. 
As for equation (2.1), we assume that b and c > are such that r > almost surely and such that 
(2.3) has a unique solution; see, for example, Karatzas and Shreve (1991, Section 5.2). 

Prom the principle of no-arbitrage in the bond market, there is a market price of risk process 
{qt} for the bond that is adapted to the filtration generated by {Wt}; see, for example, Lamberton 
and Lapeyre (1996) or Bjork (2004). Moreover, the bond market's price of risk at time t is a 
deterministic function of the short rate and of time; that is, qt = q{rt., t)- Thus, the time-t price F 
of a default-free, zero-coupon bond that pays $1 at time T is given by 



F{r,t;T) = 



nT 






rt = r 



(2.4) 



in which Q is the probability measure with Radon-Nikodym derivative with respect to P given 
by ^1^^ = exp J*q{r„s)dWs - \ llq^{r,,s)d^. It follows that {W^}, with Wt^ = + 

/o q{i^s,s) ds, is a standard Brownian motion with respect to Q. 

Prom Bjork (2004), we know that the price F of the T-bond solves the following PDE: 

Ft + b'^ir, t)Fr + ^c^r, t)Frr -rF = 0, 
F{r,T;T) = l, 

in which b^ = b— qc. Por the remainder of Section 2, we drop the notational dependence of F on 
T because T is understood; whenever the bond pays at a time other than T, then we specify that 
time. 

We can use the PDE (2.5) to obtain the dynamics of the T-bond's price F{rs, s), in which we 
think oi Tt = r as given and s G [t, T]. Indeed, 

f dF{rs, s) = (fs F{rs, s) + q^r^, s) c{rs, s) F^(rs, s))ds + c^r^, s) Fr{rs,s) dW^, 

\ F{rt,t) = F{r,t). ^ ' ' 

We use (2.6) in the next section when we develop the dynamics of a portfolio containing the 
obligation to pay the life annuity and a certain number of T-bonds. 
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2.2 Valuing via the Instantaneous Sharpe Ratio 

In this section, we first describe our method for valuing life annuities in this incomplete market. 
Then, we fully develop the valuation for a single life annuity, that is, n = 1, and we present the 
valuation for the case when n > 1. 

2.2.1 Recipe for valuation 

The market for insurance is incomplete; therefore, there is no unique pricing mechanism. 
To value contracts in this market, one must assume something about how risk is "priced." For 
example, one could use the principle of equivalent utility (see Zariphopoulou (2001) for a review) 
or consider the set of equivalent martingale measures (Blanchet-Scalliet, El Karoui, and Martellini, 
2005; Dahl and M0ller, 2006) to value the risk. We use the instantaneous Sharpe ratio because 
of its analogy with the bond market's price of risk and because of the desirable properties of the 
resulting value; see Sections 3.1 and 4. Also, in the limit as the number of contracts approaches 
infinity, the resulting value can be represented as an expectation with respect to an equivalent 
martingale measure. Because of these properties, we anticipate that our valuation methodology 
will prove useful in pricing risks in other incomplete markets; for example, see Bayraktar and 
Young (2008). 

Our method for valuing contingent claims in an incomplete market is as follows: 

A. First, define a portfolio composed of two parts: (1) the obligation to underwrite the contingent 
claim (life annuity, in this case), and (2) a self-financing sub-portfolio of T-bonds and money 
market funds to (partially) hedge the contingent claim. 

B. Find the investment strategy in the bonds so that the local variance of the total portfolio 
is a minimum. If the market were complete, then one could find an investment strategy so 
that the local variance is zero. However, in an incomplete market, there will be residual 
risk as measured by the local variance. This control of the local variance is called local risk 
minimization by Schweizer (2001a). 

C. Determine the value of the contingent claim so that the instantaneous Sharpe ratio of the 
total portfolio equals a prc-spccified value. This amounts to setting the drift of the portfolio 
equal to the short rate times the portfolio value plus the Sharpe ratio times the local standard 
deviation of the portfolio. Cochrane and Saa-Requejo (2000) and Bjork and Slinko (2006) 
apply the idea of limiting the instantaneous Sharpe ratio to restrict the possible prices of 
claims in an incomplete market, and in Section 3.1, we further discuss how our method relates 
to theirs. Also, Schweizer (2001b) derives a financial counterpart to the standard deviation 
premium principle. His was a global standard deviation premium principle, whereas ours can 
be viewed as a local standard deviation premium principle. 

2.2.2 Valuing life annuities 

Consider the time-t value of the life annuity that pays at the continuous rate of $1 per year 
until the individual dies or until time T, whichever comes first. Denote the value of this annuity 
by a = a{r,X,t), in which we explicitly recognize that the value of the annuity depends on the 
short rate r and the hazard rate A at time t. (As an aside, by writing a to represent the value of 
the annuity, we mean that the buyer of the annuity is still alive. If the individual dies before time 
T, then the value of the policy jumps to 0. Also, the value at time T of the annuity is 0.) 
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As prescribed in Step A of Section 2.2.1, suppose the insurer creates a portfolio with value 
lit at time t. The portfolio contains (1) the obligation to underwrite the life annuity, with value 
—a, including the continuous flow of $1 out of the portfolio to pay the life annuity, and (2) a self- 
financing sub-portfolio of T-bonds and money market funds with value Vt at time t to (partially) 
hedge the risk of the annuity. Thus, while the individual is alive, lit = —o-iftj^tjt) + Vt with a 
continuous outflow of money at the rate of $1 per year for t < T. Let TTt denote the number of 
T-bonds in the self-financing sub-portfolio at time t, with the remainder, namely Vt ~ T^tFiTti't)^ 
invested in a money market account earning the short rate r^. 

The annuity risk cannot be fully hedged because of the randomness inherent in the individual 
living or dying. If the individual dies at time t < T, then the value of the annuity jumps to 0; 
therefore, the value of the portfolio Ht jumps from --a^rt, Xt,t) + Vt to Vt- In other words, the 
value of the portfolio instantly changes by a{rt, Xt,t). Note that the portfolio that includes the 
value of the annuity is not self-financing, although the sub-portfolio with value Vt is. 

For t > 0, let r(t) denote the future lifetime of the individual; that is, t + T{t) is the time of 
death of the individual given that the individual is alive at time t. Alternatively, r(t) = mi{u > 
0:N{t + u) = l\N{t) = 0}. 

Proposition 2.1 When Vt = r and At = A, the drift and local variance of the portfolio lit are 
given by, respectively, 

lim ^ (E {Ut+^hAr{t)) - n* I JP-t) ) = -P^a(r, A, i) - 1 + TTt q{r, t) c(r, t) Fr{r, t) + rU, (2.7) 

and 

]im^\ar{Ut+^hAr(t))\^t)=c\r,t){TrtFr{r,,s)-ar{r,^ 

(2.8) 

Proof. We first specify the dynamics of a{rt,Xt,t) and Vt. By Ito's lemma (see, for example, 
Protter (1995)), the dynamics of the value of the life annuity are given as follows: 

da{rt, At, t) = atdt + drt + ^Qrr d[r, r]t + ax dXt + ^axx d[X, X]t 

= atdt + ar{bdt + cdWt) + ^a^r- dt + ax{Xt - X){ndt + adWt) + ]^ax\ (T^(At - Xfdt, 

(2.9) 

in which we suppress the dependence of the functions a, fe, etc. on the variables r. A, and t. Also, 
[r, r], for example, denotes the quadratic variation of r. Because the sub-portfolio with value 
Vt = 'KtF{rt,t) + {Vt — TrtF{rt,t)) is self-financing, its dynamics are given by 

dVt = TTt dF{rt, t) + rt{Vt - TrtF{rt, t))dt ={TTtqcFr + rt Vt) dt + TTt cF^ dWt, (2.10) 

in which the second equality follows from (2.6), and we again suppress the dependence of the 
functions on the underlying variables r and t. 

Prom the dynamics in equations (2.9) and (2.10), from the fact that the portfolio continually 
pays at the rate of $1 while the individual is alive, and from the jump in the portfolio value when 
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the individual dies (see the discussion preceding Proposition 2.1), it follows that the value of the 
portfolio at time t + {h A T{t)) with h> 0, namely ^t+{hAT{t)): equals 

pt+{hAT{t)) pt+{hAT{t)) nt+{hAT{t)) 

nt+(?iAT(t)) = n* - / da{rs,Xs,s) + dV^ - ds 

Jt Jt Jt 

I't+ihArit)) 

+ J a{rs,Xs,s)dNs 

l-t+{hAT{t)) pt + {hAT{t)) 

= lit - / {'D''a{rs,Xs,s) + l)ds + c{rs, s){7rsFr{rs, s) - ar{rs, Xs, s)) dWg 

Jt Jt 

l-t+{hAT(t)) i-t+{hAT{t)) 

- a{s){Xs - X)ax{rs,Xs,s)dW^ + a{rs,Xs,s){dNs - X^ds) 

nt+ihArit)) i-t+{hAT{t)) 

(2.11) 

in which is an operator defined on the set of appropriately differentiable functions on R"*" x 
(A,oo) X [0,r] by 

V^v = vt + bvr + ^(?Vrr + //(A - A)va + ^^^(A - A)^vaa " (r + X)v. (2.12) 

{A^t} denotes a counting process with stochastic parameter At at time t that indicates when the 
individual dies. Note that we adjust dNt in (2.11) so that Nt — J^^Xgds is a martingale. The 
expressions in (2.7) and (2.8) follow from (2.11). □ 

In this singie-lifc case, the process 11 is "killed" when the individual dies. If we were to con- 
sider the value a'^"^ of n conditionally independent and identically distributed lives (conditionally 
independent given the hazard rate), then N would be a counting process with stochastic parameter 
nAt until the first death such that 11 jumps by a(") - a("-i) upon that death, and N becomes a 
counting process with stochastic parameter (n — l)At, etc. We consider a^"^^ later in this section 
and continue with the single-life case now. 

As stated in Step B in Section 2.2.1, we choose the number of T-bonds nt to minimize the 
local variance of this portfolio, namely lim/j^o ji^^''^(^t+{hAT{t))\^t), a dynamic measure of risk of 
the portfolio. Note that this investment strategy is the same one advocated by Schweizer (2001a) 
in local risk minimization. We have the following corollary of Proposition 2.1. 

Corollary 2.2 The optimal investment strategy ttj that minimizes the local variance is given by 

< = ar{rt,Xt,t)/Fr{rt,t), 
and under this assignment, the drift and local variance become, respectively, 

hm i (E {U^+^hArit)) l^t) - n) = -V'^'^air, X,t)-1 + rU, (2.13) 

and 

lim ivar {Ut+^hAr{t))\^t) = ^r^t^X - Xfalir, A, t) + Xa\r, X, t), (2.14) 

in which the operator V^'^ is defined in (2.12) with b replaced by b^ = b — qc. □ 
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Now, we come to valuing this annuity via the instantaneous Sharpe ratio, as in Step C of 
Section 2.2.1. Because the minimum local variance in (2.14) is positive, the insurer is unable to 
completely hedge the risk of the life annuity. Therefore, the value should reimburse the insurer for 
its risk, say, by a constant multiple a > of the local standard deviation of the portfolio. It is this 
a that is the instantaneous Sharpe ratio. We restrict the choice of a to be bounded above by 
to avoid values that violate the principle of no arbitrage; see Section 3.1 for further discussion of 
this point. 

To determine the value a of the life annuity, we set the drift of the portfolio equal to the short 
rate times the portfolio plus a times the local standard deviation of the portfolio. Thus, from 
(2.13) and (2.14), we have that a solves the equation 



-V^ a- 1 + rn = rn + a^a^{t){X - X^al + Xa^. (2.15) 
We summarize the above discussion in the following proposition. 

Proposition 2.3 The value of the life annuity a = a{r, A, t) solves the non-linear PDE given by 

at + b'^Ur + ^c^Qrr + At(A - X)ax + ^'^^('^ ~ A)^aAA - (r + A)a + 1 

= -a^(72(A-A)2a2 +Aa2, ^2-^^) 
a(r,A,r) =0, 

in which the terminal condition arises from the assumption that the annuity only pays until time 
T. □ 

If we had been able to choose the investment strategy tt so that the local variance in (2.14) 
were identically zero (that is, if the risk were hedgeable), then the right-hand side of the PDE in 
(2.16) would be zero, and we would have a linear differential equation of the Black-Scholes type. 

To end Section 2.2.2, we present the PDE solved by the value a^") of n conditionally indepen- 
dent and identically distributed life annuity risks. Specifically, we assume that all the individuals 
are of the same age and are subject to the same hazard rate as given in (2.1); however, given that 
hazard rate, the occurrences of death are independent. As discussed in the paragraph following 
the proof of Proposition 2.1, when an individual dies, the portfolio value 11 jumps by 
By paralleling the derivation of (2.16), one obtains the following proposition. 

Proposition 2.4 The value (r. A, t) of n conditionally independent and identically 

distributed life annuity risks solves the non-linear PDE given by 

' af^ + 6^4") + )^c^af) + ii{X - A)ai") + ^^^(A - Xfa^;^l - ra^^^ - nX (a^") - a^"-^)) + n 

[a(")(r,A,r)=0 

(2.17) 

The initial value in this recursion is a^^^ = 0, and the value a defined by (2.16) is a^^\ 

Our model includes the special case for which a = 0; this corresponds to the local risk 
minimization pricing and hedging model of Schweizer (2001a). Therefore, our methodology is 
more general than that of local risk minimization. 
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Corollary 2.5 When a = 0, then a^") = na°'^ solves (2.17), in which is given by 



a^\r,\t) = F{r,t-s) 



E 



Xudu 



(2.18) 



in which F is the time-t price of a default-free zero-coupon bond that pays $1 at time s as in (2.4) 
with T replaced by s. □ 



3. Relationship with the Literature 

In Section 3.1, we show that our vahiation of a life annuity for the seller is identical to the 
upper good deal bound of Cochrane and Saa-Requejo (2000), generalized by Bjork and Slinko 
(2006) to the case of a jump diffusion. The lower good deal bound is obtained as the buyer's 
valuation of the life annuity. Then, in Section 3.2, we compare our valuation with that given by 
indifference pricing via expected utility, as in Zariphopoulou (2001), for example. 

3.1 Good Deal Bounds 

Our work is this section is motivated by similar results of Bayraktar and Young (2008); sec 
the remark at the end of their Section 2.5. It is straightforward to show that one can write (2.16) 
as ^ ^ 

at + b'^ar + -^c^arr + 7:0'^(A - A)^aAA + 1 



+ « .2™^? ^ [{ii + Sa){X-X)ax-{r + X{l+^))a]=0, 
a{r,X,T) = 0. 

Note that the optimal values for S and 7 are given by, respectively, 

aa{X — X)ax 



(3.1) 



S* = 



^a^{X-X)^a;i + Xa^' 



and 
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aa 



^a^iX-X^ai + Xa^ 



(3.2) 



(3.3) 



Because we restrict a to lie between and \/A inclusively, it follows that 7* > — 1 because A > A. 
Therefore, the restriction that 7 > — 1 in (3.1) is automatically satisfied at the maximum but we 
need it for the representation of a in (3.4) below. 

Prom (3.1), it follows that a can be represented as 



a(r, A, t) = sup E 

{{5a,7j: <5i+Aa7|<a, 7.>-l, t<s<T} 



r.X.t 



e 



J'^°(''u+Au(l+7u)) du 



ds 



sup E 

{{5s,ls}- 51+Kli<cc, 7s>-l, t<s<T} 



r,X,t 



(3.4) 



Here, {rt} and {At} follow the processes dr^ = b^du-\-cdWu and dXu = {iJ,-\-Sa){Xy, — X)du-\-a{Xu — 
X)dW^, respectively, with = + Jq q{rs,s)ds and = — 5^ ds. The processes {Wt} 
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and {VF/'} are independent standard Brownian motions with respect to the filtered probabihty 
space {0,,J^, (J^t)t>o,P), in which ^1 = Lt, with Lt given by 



dLt = Lt_ {-q{rt,t) dWt + 6t dW^^ + -ft{dNt - A* dt)) , 
Lt = l. 



(3.5) 



E denotes expectation with respect to P. The expression for a in (3.4) is identical to the upper 
good deal bound given by Bjork and Slinko (2006; Theorem 2.2). 

As mentioned following the expression for 7* in (3.3), we have 7* > —1 because a < VA- 
This fact ensures that the expression in (3.4) is no greater than E*"'* [ F{r, t; s) ds] , which is 
required by no arbitrage. Also, note that 7* > — 1 is required in order for P defined via (3.5) to 
be nonnegative. 

Insurance actuaries are familiar with expressions similar to those in the expectation of (3.4) 
because in valuing life annuities, they eflfectively discount for interest and mortality often modifying 
the hazard rate to account for longevity risk (Wang, 1995). To gain further understanding of (3.4), 
we can write it as follows: 



a{r,\,t) 



F{r, t; s) sup E 

{{5,,7,}: 52+Aa7?<a, ls>-l, t<s<T} 



ds. 



(3.6) 



in which we slightly abuse notation by writing E here. 

Because As(l + 7s) > with probability 1, we can think of it as a modified hazard rate in 
(3.6). Note that if a = 0, then 7* = and S* = 0, and the expectation in the integral of (3.6) 
reduces to y?"°(A,t;s) := E'*''* [ex.p{— Xu du)], the physical probability that a person survives 
from time t to s, and the corresponding value of a equals o"'' as in (2.18). 

One can repeat the derivation in Section 2.2 to show that the value a!^ for the buyer of a 
life annuity solves (2.16) with a replaced by —a. Equivalently, solves (3.1) with max replaced 
by min. Therefore, the buyer's value of the annuity is identical to the lower good deal bound of 
Bjork and Slinko (2006; Theorem 2.2). It follows that we have the following bid-ask interval for 
life annuity values under our methodology, which we have just shown is equal to the one under 
good deal bounds: 



inf E 

{5a ,7s} 



r,A,t 



■ f rftt -, 7 

i{T{t)>.s}ds 



, sup E""'^'* 

{<5s;7s} 



^{T{t)>s}ds 



(3.7) 

in which we take the sup, respectively in/, over {{5s, 7s} : ^1 + As 7s < a, 7s > —1, t < s < T}. 
Therefore, we have provided an alternative derivation of the good deal bounds for our setting as 
the result of a local risk minimization investment strategy combined with a (local) risk loading 
expressed as a multiple of the local standard deviation of the portfolio. 

The interval in (3.7) is a subinterval of the interval of no-arbitrage prices 



inf E*^ 

QeM 



1: 



i: 



r-u du 



l{T(t)>s} ds 





, sup E*^ 


[i 




QeM 





r(t)>4 



ds 



(3.8) 



in which M. is the set of equivalent martingale measures. Since this no-arbitrage pricing interval is 
too wide and practically useless, it was Cochrane and Saa-Requejo (2000) 's idea to find a tighter 
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subinterval. From the representation of the seller's price in (3.4), one can see that our valuation 
method results in a coherent risk measure (Artzner et al., 1999). 



3.2 Indifference Pricing via Expected Utility 



Ludkovski and Young (2008) apply indifference pricing via expected utility to value pure 
endowments and life annuities in the presence stochastic hazard rates. They use exponential 
utility so that the value of the insurance contracts are independent of the wealth of the seller. 
With absolute risk aversion r/ > 0, one can follow the derivation of equation (2.41) in Ludkovski 



and Young (2008) to show that, in the setting of this paper, the indifference price a^^ 
solves the following PDE: 



^ip 



{r,X,t) 



1 



1 



2„IP 



2F 



-Tja 



xx-yr + A)a^^ + 1 
IP • 



(3.9) 



(r,A,r) = 0. 



Here F = F{r,t;T). Note that when rj = 0, then a^^ = a"°, as defined in (2.18), which is also 
the price under the local risk minimization method of Schweizer (2001a). The choice of a in (2.16) 
reflects the seller's risk aversion, just as the choice of rj in (3.9) reflects the seller's risk aversion in 
the setting of expected utility. 

Parallel to equation (3.1), one can rewrite (3.9) as follows: 



ar + b<^ai^ + Ic'ai^ + '-a^X - Xfai^, + 1 



+ max 

5, 7 



F \F 
-ir'^' + — (7-(l + 7)ln(l + 7)) 



(3.10) 



+ (/X + Sa){X - X)ai^ - (r + A(l + j))a 



JP 



0, 



U''^(r,A,r) = o. 



From (3.10), it follows that a^^ can be represented as 



a^^(r,A,t) 



sup E 

{<5s,7a} 



r.X.t 



-f +A.(7.-(l + 7.)ln(l + 7«)) 



ds 



(3.11) 

As in (3.4), {rt} and {At} follow the processes rir,j = iPdu + cdWu and dXu = {jJL + 5a){Xu — 
X)du + a{Xu - X)dW^, respectively, with Wu = Wu + q{rs, s)ds and = - 5„ ds. Note 
that (3.11) is similar to the representation of a in (3.4) with the constraint on the controls {(5s,7s} 
replaced by the penalty term {F{rs, s;T)/r]) {-6^/2 + Xs{'js - (1 + 7s)ln(l + 7s))). 

What follows is a formal discussion to show how the PDE in (3.9) is further related to the one 
in (2.16). For rja^^ /F small enough, one can "approximate" this PDE by expanding e"''" and 
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ignoring terms of the third power and higher. When one does this, the resulting PDE is given by 

' At + }fiAr + ^c^Arr + /x(A - X)A^ + ^a\\ - A)2Aaa - (r + A)^ + 1 
' =-^{-H^-^?Al + XA^), (3-12) 
.^(r,A,r) = 0. 

Compare this PDE with the one in (2.16). Note that they are quite similar with a in (2.16) 
replaced by r?/(2F) in (3.12) and with the removal of the square root in (2.16). 

Milevsky, Promislow, and Young (2005) show that valuing insurance risks via the instantaneous 
Sharpe ratio is closely related to the standard deviation principle, which is used by actuaries to 
price insurance (Bowers et al., 1986). Similarly, Pratt (1964) shows that indifference pricing is 
closely related to the variance principle, under which the risk load is the square of the risk load 
under the standard deviation pricing principle. Pratt's result is consistent with (3.12). Thus, the 
valuation method that we propose in this paper is related to indifference pricing in a way parallel 
to how the standard deviation principle is related to the variance principle, both of which are used 
in insurance pricing. 

4. Properties of a^"^ the Value of n Life Annuities 

To demonstrate important properties of a^'^\ we rely on a comparison principle (Walter, 1970), 
and we present that in Appendix A. In Section 4.1, we list some qualitative properties of a*^"^ . Then, 
in Section 4.2, we prove the limiting result for ^o^"-*- 

Assumption 4.1 Henceforth, we assume that the volatility on the short rate c satisfies the growth 
condition in the hypothesis of Theorem A.l and that the drifts and fj, satisfy the growth 
conditions in the hypothesis of Lemma A. 2. For later purposes (specifically for Property 3 in 
Proposition 4.1 below), we also assume that fix satisfies the growth condition |ma|(A — A) + < 
K (^1 + (ln(A — A))^^ . Also, we assume that < a < ^/X, as stated in Section 3.1. 

4-1 Qualitative Properties of a^"^^ 

In this section, we list and discuss properties of the risk-adjusted value a^"^ of n life annuity 
contracts. We state the following proposition without proof because its proof is similar to those in 
Milevsky, Promislow, and Young (2005). 

Proposition 4.1 Under Assumption 4.1, a(") = a(") (r. A, t) satisfies the following properties on 
G = R+ X (A, oo) X [0, T] forn>l: 

(1) No arbitrage: < o^") < n F {r, t; s) ds . 

(2) Increasing in n: a*-"^ > a^"'~^\ 

(3) Decreasing in A; a^"'' < 0. 

(4) Increasing in a: Suppose < ai < 02 < and let a'-"'''"' be the solution of (2.17) with 
a = ai, for i = l,2 andn>0. Then, a(")'«i < a^'^^ . 

(5) Lower bound: na°''^ < a*^")'", in which is defined in (2.18). 

(6) Decreasing in ji: Suppose ni < fi2, and let a^"^''** denote the solution of (2.17) with ji = Hi, 
fori = 1,2. Then, a(")'''i > a^")'''^ 
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(7) Increasing in a if convex in X: Suppose < < a2{t) on [0, T], and let a^"-')''^' denote the 
solution of (2.17) with a = ai, for i = 1,2. If a^^^''^^ > for all n, or if a^^^''^^ > for all n, 
then aW''^i < a^''^^ 

(8) Subadditive: For m, n nonnegative integers, a(™+") < a(™) + a^""^ . 

(9) Decreasing value per risk: ^a^"'^ decreases with respect to n>\. 

(10) Scaling: The value of n annuities that pay at a rate of k > is ka^"'\ 

Remarks: 

4.1 Because the payoff under a life annuity is nonnegative, we expect its value to be nonnegative. 
Also, if we hypothesize that the individuals will not die, then we get the upper bound given 
in Property 1. In Section 4.2, we sharpen the bounds given in Property 1 considerably. 

4.2 One uses Property 2 to prove Property 3. However, Property 2 is interesting in its own right 
because it confirms our intuition that the marginal value of adding an additional policyholder, 

namely a*^"^ — a^"~^\ is nonnegative. 

4.3 Property 3 makes sense because we expect the value of life annuities to decrease as the hazard 
rate increases, that is, as individuals are more likely to die. 

4.4 Property 4 states that as the parameter a increases, the risk-adjusted value a^"'^'°' increases. 
This result justifies the use of the phrase risk parameter when referring to a. Also, by referring 
to the expression in (2.18), we see that Property 5 follows from Property 4. Therefore, the 
lower bound of i-a^")'" (as a approaches zero) is the same as the lower bound of a" = a^^^'", 
namely, . We call the difference > the risk charge per annuity for a portfolio 
of n life annuities. 

Because > 0, the strategy we propose results in a net profit on average, but 

this margin is necessary to compensate the issuer for the undiversifiable mortality risk, as in 
indifference pricing (Section 3.2). 

4.5 In words, Property 6 tells us that as the drift of the hazard rate increases, then the value of the 
life annuities decreases. This occurs for essentially the same reason that the value decreases 
with the hazard rate; see Property 3. 

4.6 Subadditivity is a reasonable property because if it did not hold, then buyers of annuities 
could buy their annuities separately and thereby save money. Subadditivity follows as a 
corollary from Property 9, which is interesting in itself, namely that the average value per 
risk decreases as the number of risks increases. In other words, the risk charge per annuity, 
namely > 0, decreases as n increases. In the next section, we will determine the 
limiting value of the risk charge per annuity as n goes to infinity, and from this result, we 
will subdivide the risk charge into the portion due to the finite size of the portfolio and the 
portion due to randomness of the hazard rate. 

4.7 The same ten properties hold for the value of n pure endowment risks (given in Appendix B), 
as shown in Milevsky, Promislow, and Young (2005), with Property 1 appropriately modified. 

These properties-which also hold in simpler, static settings (that is, for typical actuarial 
premium principles)-help to demonstrate that our valuation mechanism is reasonable. What we 
gain in going to our more complicated setting is a dynamic pricing mechanism that compensates 
the issuer of the annuity for the undiversifiable mortality risk and shows the issuer how to hedge 
using a local risk minimization strategy. 
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4-2 Limiting Result for ^o^"^ 

In this section, we determine and interpret the hmit hm^^cx) In the following two 

lemmas, we apply Theorem A.l and Lemma A. 2 to show that on G = R"*" x (A, oo) x [0,T], we 
have n F{r,t; s) (]{X,t; s) ds < a(")(r,A,t) < F{r,t; s) (p^'''> {X,t; s) ds, in which /3 and i^^") 
are defined in (B.4) and (B.2), respectively, in Appendix B. Thus, if we divide by n, we have 
that lim„^oo ^«*'"''(^! ^) = /* -^('"' ^) /3(A, t; s) ds. See Appendix C for the proofs of the two 
following lemmas. 

Lemma 4.2 Under Assumption 4.1, we have 

a^''\r,X,t) < j F{r,t;s)ip^"'\\,t;s)ds, (4.1) 

on G for n> 0, in which t/j^") is defined in (B.2). 

Prom Lemma 4.2, we learn that if one were to treat (/j*^"^ as a modified expected number of 
survivors in a traditional actuarial computation of an annuity value, that is, the right-hand side of 
(4.1), then we overvalue the annuity. Inequality (4.1) makes sense because a^'^'> takes into account 
the reduction in overall risk when valuing an annuity (left-hand side of (4.1)) versus a series of 
essentially independent infinitesimal pure endowments (right-hand side of (4.1)). 

Lemma 4.3 Under Assumption 4.1, we have 

fT 

nj F{r,t;s)P{X,t;s)ds<a^''\r,X,t), (4.2) 

on G for n>0, in which f3 is defined in (B.4), or equivalently, in (B.6). 

Lemma 4.3 implies that the value a*^""-* of a life annuity is greater than n times the value of a life 
annuity computed by using the limiting value of -(^(") , namely (3. In particular, a(") is greater than 
n times the traditional actuarial value of a life annuity, namely ff F{r, t; s) <^"°(A, t; s) ds because 
/3 is greater than or equal to the physical probability of survival ip"^. Recall from Appendix B that 
the value of n conditionally independent and identically distributed pure endowments is given by 
F(fi^"-\ 

We combine these two lemmas with a result from Milevsky, Promislow, and Young (2005) to 
obtain the main result of this paper. 

Theorem 4.4 For (r. A, t) G G, 

1 r'^ 

lim -a(")(r, A, t) = / Fir,t; s) P{X,t; s) ds. (4.3) 

Proof. Prom Milevsky, Promislow, and Young (2005), we know that the sequence {^f^"'^) de- 
creases uniformly to (3 on (A, oo) x [0, s] for all < s < T. Also, recall from Property 9 in Propo- 
sition 4.1 that ^ot^"^ decreases with respect to n. Because F(pds < (pds < - — < oo, 

it follows from the Lebesgue (dominated) convergence theorem (Royden, 1988, page 267) and 
Lemmas 4.2 and 4.3 that 
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/ Fpds< Urn -o(") < [ F lim - = [ F/3ds, (4.4) 

Jt n Jt it 

which implies (4.3). □ 

Prom the discussions accompanying (2.4) and (B.6), we can represent the Umit in (4.3) as an 
expectation. 

Corollary 4.5 For (r, A, t) G G, 



hm -a(")(r, A, t) =E'^'^'* / e'T^ 

n—^oo n J I 



(4.5) 



in which {rt} and {A*} follow the processes dru = h^du + cdWu and dXu = {ix — aa){\u — X)du + 
^■(Au — K)dW^, respectively, with = + q(rs,s)ds and = + au. The processes 
{Wt} and {W^} are independent standard Brownian motions with respect to the filtered probability 
space {n, {J't)t>o, P), in which § = exp (- £ q{r,, s) dW, - | q''{r,,s) ds) e-^'^'-Wt. 
E denotes expectation with respect to P. 
Remarks: 

4.8 The expression in (4.5) is identical to that of the value of a coupon bond that pays coupons 
continuously at the rate of $1 under the discount rate of rj + Aj, instead of simply rj as for a 
bond that cannot default. Therefore, one can think of At as the default rate. For this reason, 
we anticipate that our methodology will prove useful in valuing credit risk derivatives. 

4.9 The expression in (4.5) is an expectation with respect to a measure that is equivalent to the 
physical measure P. Therefore, in the limit, we obtain an arbitrage-free value for the life 
annuity. Thereby, one can think of oTir valuation method as generalizing the one of Blanchet- 
Scalliet, El Karoui, and Martellini (2005). 

4.10 If the hazard rate is deterministic, then in the limit, there is no mortality risk. Specifically, if 
(7 = 0, then lim„_>oo ■^ot^"^ = on G, in which a"" is the net premium given in (2.18). In 
other words, if the hazard rate is deterministic, then as the number of contracts approaches 
infinity, the value of a life annuity collapses to the discounted expected payout using the 
physical probability measure to value the mortality risk, regardless of the target value of the 
Sharpe ratio. Alternatively, one can see that in the limit the average cash fiow is certain, 
hence, the value becomes that of a "certain" (or non-life) annuity with a rate of discount 
modified to account for the rate of dying. 

4.11 On the other hand, if the hazard rate is truly stochastic, then in the limit, mortality risk 
remains. Specifically, if a is uniformly bounded below by k > 0, then lim^^co n^^"^^ ^ ^"'^ on 
G, with equality only when t = T. 

4.12 We are now prepared to subdivide the risk charge per annuity into the portion due to the finite 
portfolio and the portion due to the stochastic hazard rate. Prom the proof of Theorem 4.4, 
we know that ^a'^'^^r, X,t) > F(r,t;s) f3{X,t;s)ds, with equality in the limit. Therefore, 
define ^o*^"-* — F (3 ds as the risk charge (per risk) for holding a finite portfolio of n annuities, 
and define F /3ds — as the risk charge for the stochastic hazard rate. Thus, we have 

ia<») - a"' = (la<-.) - r^-S*) + ( I^F^is - a"> \ , (4.6) 
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in which the risk charge for the stochastic hazard rate is zero if cr = and is positive otherwise. 

4.13 The representation of the hmiting value of ia^"^) in (4.3) is quite similar to a net single 
premium as defined in courses on Life Contingencies (Bowers et al., 1986). Indeed, F is the 
monetary discount function, and f3 plays the role of the survival probability. Additionally, we 
can express the limit F(r, t; s) t; s) ds as the solution p of the following PDE: 

Pt + b^pr + ^c^Prr + (m " ^a) (A - X)p^ + (A - Xfpxx - (r + X)p + 1 = 0, 
p{r,X,T)=0. 

Note that a is analogous to the bond market's price of risk q in (2.4); therefore, we refer to a 
as the (annuity) market's price of mortality risk. 

4.14 The proof of Theorem 4.4 suggests that we can use other reasonable models for bond pricing. 
Additionally, one could use a different sort of default-free bond, such as a consol bond (or 

perpetuity). In fact, in work not shown here, the authors obtained the same PDEs as in 
(2.16) and (2.17) when replacing the T-bonds with consol bonds. The latter may be more 
appropriate in the setting of whole life annuities. Alternatively, an insurer could continually 
roll money into longer-term bonds as the shorter-term bonds mature. 

4.15 By comparing Corollary 4.5 with the expression in (3.4) and the discussion afterwards, we see 
that as n approaches infinity, the controls in the expectation in (3.4) are further restricted so 
that 7s = with probability 1, that is, 5s = —a with probability 1. 

4.16 Only in the limit does our valuation method result in a linear pricing rule, as evidenced by the 
linear PDE in (4.7) as compared with the non-linear PDE in (2.16). Pricing via an equivalent 
martingale measure (EMM) results in a linear pricing rule; therefore, our general valuation 
method is not equivalent to choosing an EMM. 

4.17 When valuing life insurance using the method of this paper, Young (2007) proves results 
parallel to those in Sections 4.1 and 4.2. In particular, it is encouraging to learn that the 
limiting results are robust to the type of insurance contract (life insurance versus life annuities). 
Also, in that paper, there is a numerical algorithm for valuing life insurance. An interested 
reader can easily modify that numerical algorithm to value life annuities. 

5. Conclusion 

A number of recent research papers focused attention on the valuation of longevity risk from 
a variety of perspectives. In this paper, we returned to basics and showed how the uncertainty in 
the hedging portfolio translates into a non-zero standard deviation per policy. We draw upon the 
financial analogy of the Sharpe ratio to develop a methodology for valuing the non-diversifiable 
component of aggregate mortality risk. Our main qualitative insight is that similar to the financial 
economic approach to analyzing stock market risk, the uncertainty embedded within mortality- 
contingent claims can be decomposed into idiosyncractic (diversifiable) and non-diversifiable com- 
ponents. 

We developed a theoretical foundation for valuing mortality risk by assuming that the risk is 
"priced" via the instantaneous Sharpe ratio. Because the market for life annuities is incomplete, 
one cannot assert that there is a unique price. However, we believe that the price that our method 
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produces is a valid one because of the many desirable properties that it satisfies and because it is 
the upper good deal bound of Cochrane and Saa-Requejo (2000), generalized to the jump diffusion 
setting by Bjork and Slinko (2006). Therefore, we have the additional contribution of providing 
another derivation of the good deal bounds for the setting of this paper, one that includes a hedging 
strategy along with a actuarially-inspired risk loading in the form of the local standard deviation. 

We also studied properties of the value of n conditionally independent and identically dis- 
tributed life annuities. The value is subadditive with respect to n, and the risk charge per person 
decreases as n increases. We proved that if the hazard rate is deterministic, then the risk charge 
per person goes to zero as n goes to infinity. Moreover, we proved that if the hazard rate is 
stochastic, then the risk charge person is positive as n goes to infinity, which reflects the fact that 
the mortality risk is not diversifiable in this case. Additionally, in Remark 4.12, we decomposed 
the per-risk risk charge into the finite portfolio and stochastic mortality risk charges. 

Appendix A. Comparison Principle 

In this appendix, we present a comparison principle from Walter (1970, Section 28) on which 
we rely extensively in the proofs of the properties of a^"^ given in Section 4. We begin by stating 
a relevant one-sided Lipschitz condition along with growth conditions. Suppose a function g = 
g{r,X,t,v,p,q) satisfies the following one-sided Lipschitz condition: For v > w, 

g{r,X,t,v,p,q)-g{r,X,t,w,p' ,q') < fi{r, X,t){v - w) + f2ir, \,t)\p- p'\ + fair, X,t)\q - q'\, (A.l) 
with growth conditions on /i, /2, and fs given by 

0<fi{r,X,t) <i^(l + (lnr)2 + (ln(A-A))2), (A.2a) 

0</2(r,A,t) <i^r(l + |lnr| + |ln(A-A)|), (A.26) 

and 

0</3(r,A,t) <K(A-A)(l + |lnr| + |ln(A-A)|), (A.2c) 

for some constant K > 0, and for all {r,X,t) G R"*" x (A, oo) x [0,T]. To prove properties of a^"-\ 
we use the following comparison principle, which we obtain from Walter (1970, Section 28). For 
the proof of a similar result, see Milevsky, Promislow, and Young (2005). 

Theorem A.l Let G = R"*" x (A, oo) x [0,T], and denote by Q the collection of functions on G 
that are twice- differentiable in their first two variables and once-differentiable in their third. Define 
a differential operator C on Q by 

1 1 

ji:V = Vt + -C^{r, t)Vrr + 2^^(0(^ " Kfvxx + 9{r, A, t, V,Vr,Vx), (A.3) 

in which g satisfies (A.l) and (A.2). Suppose v,w G G are such that there exists a constant K >0 
withv< e^{(i»0'+(MA-A))^} andw> -e^lO'iO'+CMA-A))^} for large (lnr)2 + (ln(A- A))^. Also, 

suppose that there exists a constant K > such that c{r,t) < Kr for all r > and < t < T. 
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Then, if (a) Cv > Cw on G, and if (b) v{r, A, T) < w{r, A, T) for allr > and A > A, then v < w 
on G. 

As a lemma for results to follow, we show that the differential operator associated with our 
problem satisfies the hypotheses of Theorem A.l. 

Lemma A. 2 Define Qn, for n> 1, by 

g^ir, A, t, v,p, q) = b^{r, t)p + /x(A, t)(A - X)q - rv + n - nX (v - a^""^M 

I i ^^-'^ 

+ a^a^{t){X - A)2g2 + nX{v- a("-i)) , 

in which a^^~^^ solves (2.17) with n replaced by n — 1. Then, g-n satisfies the one-sided Lipschitz 
condition (A.l) onG. Furthermore, if\h^{r,t)\ <Kr(l + |lnr|) and\ii{X,t)\ < iir(l + | ln(A — A)|), 
then (A.2) holds. 

Proof. Suppose v > then 



gn{r,X,t,v,p,ci) - gn{r,X,t,w,p' ,q') 

{r,t){p - p) + fi{X,t){X - X){q - q) - (r + nX) {v - w) 



+ ay a2(t)(A - A)2g2 +nX{v- a("-i)) 

- a^a^{t){X - A)2(g')2 + nX{w- a("-i)) 

< \b'^{r,t)\\p-p'\ + {\f^iX,t)\ + aa{t)){X-X)\q-q'\- (r + nX-aV^^ {v - w) 

< Ib'^ir, t)\\p- p'\ + MX,t)\ + aa{t)) (A - A)|g - q'\. 

In the first inequality, we use the fact that A> B, then VC^ + A'^ - VC^ + B'^ < A - B. For 
the second inequality, recall that < a < \/A. Thus, (A.l) holds with /i(r, A,t) = 0, f2{r,X,t) = 
\bQ {r,t)\, and f sir, X,t) = |/Lt(A,t)|+acr(t). Note that /a satisfies (A.2b) if |6'3(r, i)| < Kr(l+|lnr|), 
and /a satisfies (A.2c) if |/x(A, t)\<K{l+\ ln(A - A)|). □ 

Appendix B. Limit Result for Valuing n Pure Endowments 

In this appendix, we review some of the results of Milevsky, Promislow, and Young (2005). Via 
the instantaneous Sharpe ratio, they value a pure endowment that pays 1 at time s if an individual 
is alive at that time. They also study properties of the time-t value P^") = P(") (r, A, t; s) of n 
conditionally independent and identically distributed pure endowment risks with <t < s. 

By paralleling the derivation of (2.17), one can show that P^") solves the non-linear PDE 
given by 



' P/") + 6^P,(") + ^c^pW + /.(A - A)Pf ) + ^a^iX - A)2p}^) - rP^ 



- nX (p(") - p('^-i)) = -a^a2(A- A)2 (p^"))\nA (pW - p(n-i))' 
^pW(r,A,s;s) = n, 



(B.l) 
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with = 0. 

We can multiplicatively separate the variables r and A in P^"). Indeed, P^"'\r, X,t; s) 
F{r,t;s) (p^'^\X,t;s), in which (p^"^^ solves the recursion 



= -a^a^X - A)2 ((^J^)) ' + nA - 
^^^^\X,s;s) = n, 



(B.2) 



with (/?(°) = 0. We can interpret (p^"'^ as a risk-adjusted expected number of survivors to time s 
from the n individuals alive at time t. Indeed, Milevsky, Promislow, and Young (2005) show that 



nE 



< < ne"(^""VA)(*"*\ and (^^"^ < 0. The main result of that paper is 



ip''"' s ne \- "V-/V" ana (p^^ 
given in the following theorem: 
Theorem B.l Under Assumption 4.1, we have 



lim -P(")(r,A,f;s) = F{r,t; s) (3{X,t; s), 



n—>oo n 

uniformly on R"*" x (A, oo) x [0, s] . Here (3 is the solution of 

(3t + {l^- a<j){X - X)(5x + ^t^'(A - Xfpxx - A/3 = 0, 
JiX,s;s) = l. 

More precisely, Milevsky, Promislow, and Young (2005) show that 



(B.3) 



(B.4) 



1 



PiX,t;s)<-p^^\X,t;s), 



n 



(B.5) 



with equality in the limit. Note that solves a linear PDE; therefore, one can represent it as an 
expectation via the Feynman-Kac Theorem (Karatzas and Shreve, 1991): 



/3{X,t;s) = E 



Xudu 



(B.6) 



in which {Aj} follows the process dA„ = {i^ — aa){Xu — X)du + a{Xu — X)dW^ , with = W^ + au. 
The process {W'^^j is a standard Brownian motion with respect to the filtered probability space 
(.7^t)t>o, P), in which = exp (^—aW^ — ^oi^t). E denotes expectation with respect to 

P. 

C. Proofs of Lemmas 4.2 and 4.3 

In this appendix, we supply the proofs of Lemmas 4.2 and 4.3. 

C.l Proof of Lemma 4-2 

We proceed by induction. The inequality is clearly true for n = because both sides equal 
zero in that case. Assume that inequality (4.1) holds for n — 1 in place of n and show that it holds 
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for n. Define a differential operator £ on ^ by (A. 3) with g = in (A.4). Because a^") solves 
(2.17), it follows that Ca^"-^ = 0, and 



C J F{r, t; s)^(") (A, t; s)ds = -F{r, t; t)^^^^ (A, i; t) + J [Ftip^^^ + F^^^^ j ds 
+ ^(T^iX-X)J^ Fip'-^^^ds-r Fif^^-^ds-nxlJ^ - a^""^) j + ; 



(C.l) 



+ a 



\ 



C72(A-A)2^^ Fifi'^^^Us^ ^"^(/ F<f(^)ds-d^-^')'^ . 



Because F{r, t;t) = 1 and (p'^'^^X, t; t) = n, the first term on the right-hand side of (C.l) cancels the 
n at the end of the third line. Also, we use the fact that F solves Ft+b^{r, t)Fr + ^c^{r, t)Frr-rF = 
and that (p^"'^ solves (B.2) to simphfy (C.l) and thus obtain 

cj^ F{r,t;s)^^''\X,t;s)ds = -nx(^J^ F^^'-^^ds - a^'^"^)^ 



-aj^ F^a^X-X)^ (ypj^^) '+ nA (</?(")- ^(""i))' 



ds 



(C.2) 



+ Q!^ (72(A- A)2 ^ Fip^^Usj F^Wc/s - a("-i) 

Now, C (/f F^p^-^Ms^ < if 



pj'^ds +nX 



(C.3) 



C72(A-A)2 (^l^F^\ 

<^ F^a^X-Xy ((/?^"^)VnA {ip(^) - ifi(^-^) fds, 

where we use the assumption that < a < ^/X, which in turn implies that aVnX < nX for n > 1 
and A > A. 

Recall that the Minkowski inequality tells us that if < p < 1, then + H^Hp < ||/ + dWp, 
in which = (/ |/|^dz/)^^^ for some measure v, (Hewitt and Stromberg, 1965, page 192). Let 

p = 1/2, / = (72(A - A)2 (^ip^^^^ , g = nX ((^W - (p("-^))^ and du = F{r,t;s)ds. Then, the 
Minkowski inequality implies 
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> ^^^(A - Xf + nX I j^^ F (^(") - ^('^-^)) 

Thus, inequality (C.3) holds if 



(C.4) 



C72(A-A)2^j^ Fif'^^'^ds^ +nx(^J^ - a("-i) j 



(C.5) 



< 



If wc define Bx = (j{t){X - A) // F ip^-^> ds, A = V^J^ F ^p^""^ ds, C = ds, and 

D = \/nXa^'^~^\ then (C.5) is equivalent to 

^Bl + {A- -{C-D)< ^Bl + {A- Cy. (C.6) 

Prom Milevsky, Promislow, and Young (2005), we know that A > C. Prom the induction as- 
sumption, we have that C > D. Thus, the left-hand side of (C.6) is positive, and demonstrating 
inequality (C.6) is a straightforward matter of squaring both sides and simplifying. 

Thus, we have that c(^J^ F ds^ < = Ca'^''\ Also, both F ip^"^"! ds and a^") equal 

when t = T. It follows from Theorem A.l and Lemma A.2 that a^") < F ds on G. □ 

C.2 Proof of Lemma 4-3 

We proceed by induction. The inequality is clearly true for n = because both sides equal 
zero in that case. Assume that inequality (4.2) holds for n — 1 in place of n and show that it holds 
for n. Define a differential operator £ on ^ by (A. 3) with g = Qn in (A.4). Because a^"'^ solves 
(2.17), it follows that >Ca(") = 0, and 



C nF{r,t;s)P{X,t;s)ds = -nF{r,t;t)P{X,t;t)+ n{FtP + FPt)ds 
+ nFrpds + ^c^ nFrrl^ds + n{X- X) nFpxds 

+ \(^'^{X-X) riFpxxds -r nFpds " ^ nFpds - a^""^) 



-I- n 



+ a 



a^{x-xy 



nFPxds^ + nX nFpds - a("-i) j 



(C.7) 



22 



Because F{r, t]t) = 1 and /3(A, t; t) = 1, the first term on the right-hand side of (C.7) cancels the n 
at the end of the third hne. Also, we use the fact that F solves Ft+b^{r, t)Fr + ^c'^{r, t)Frr — rF = 
and that /3 solves (B.4) to simplify (C.7) and thereby obtain 

cj^ nF{r,t;s)(3{X,t]s)ds = nxla^''-^'>-{n-l)J^ F^dsj 



+ an 



\ 



ct2(A-A)2 / F[5xds\ +-\ nF(3ds-a(^-^) 



n 



(C.8) 



+ ana{X-X) FPxds > = Ca^"\ 



in which the inequality follows from the induction assumption. Also, both F (5 ds and a*^"^ equal 
when t = T. It follows from Theorem A.l and Lemma A. 2 that n F P ds < a^"^^ on G. □ 
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